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ABSTRACT 

Context. Transiting highly inflated giant planets offer the possibility of characterizing their atmospheres. A fraction of the starlight 
passes through the high-altitude layers of the planetary atmosphere during transit. The resulting absorption is expected to be wave¬ 
length dependent for cloud-free atmospheres with an amplitude of up to 10“^ of the stellar flux, while a high-altitude cloud deck 
would cause a gray opacity. 

Aims. We observed the Satum-mass and Jupiter-sized exoplanet HAT-P-19b to refine its transit parameters and ephemeris as well as 
to shed first light on its transmission spectrum. We monitored the host star over one year to quantify its flux variahility and to correct 
the transmission spectrum for a slope caused by starspots. 

Methods. A transit of HAT-P-19b was observed spectroscopically with OSIRIS at the Gran Telescopio Canarias in January 2012. 
The spectra of the target and the comparison star covered the wavelength range from 5600 to 7600 A. One high-precision differential 
light curve was created by integrating the entire spectral flux. This white-light curve was used to derive absolute transit parameters. 
Furthermore, a set of light curves over wavelength was formed by a flux integration in 41 wavelength channels of 50 A width. We 
analyzed these spectral light curves for chromatic variations of transit depth. 

Results. The transit fit of the combined white-light curve yields a refined value of the planet-to-star radius ratio of 0.1390 ± 0.0012 
and an inclination of 88.89 ± 0.32 degrees. After a re-analysis of published data, we refine the orbital period to 4.0087844 ± 0.0000015 
days. We obtain a flat transmission spectrum without significant additional absorption at any wavelength or any slope. However, our 
accuracy is not sufficient to significantly rule out the presence of a pressure-hroadened sodium feature. Our photometric monitoring 
campaign allowed for an estimate of the stellar rotation period of 35.5 ± 2.5 days and an improved age estimate of 5.5^['3 Gyr by 
gyrochronology. The calculated correction of the transit depth for unocculted spots on the visible hemisphere was found to be well 
within the derived 1 cr uncertainty of the white-light curve and the spectral data points of the transmission spectrum. 

Key words, giant planet formation - /c-mechanism - stability of gas spheres 


1. Introduction 


Transiting extrasolar planets offer the opportunity of character¬ 
izing their atmospheres. During a planetary transit, parts of the 
starlight pass through an atmospheric ring at the planet’s ter¬ 
minator region, picking up the signature of the planetary at¬ 
mosphere. These spectral signatures are measurable with the 
current generation of instruments for some favorable plane¬ 
tary systems. The required specU'oscopic precision necessary 
for the detection of these features depends on the planet-to- 
star radius ratio, the scale height of the planetary atmosphere, 
and the number of photons that can be collected from the host 
star during transit. First attempts to detect additional absorp- 
tion during transit were already able to provide upper limit s 
(iRauer et alj|2000l: iBundv & Marcvl[2000l : iMoutou et al]|2001ll . 

The sensitivity neccessary for a dete ction of an atmospheric sig - 
nal was achieved for the first time bv ICharbonneau eti 


* Based on observations made with the Gran Telescopio Canarias 
(GTC), installed in the Spanish Observatorio del Roque de los Mucha- 
chos of the Instituto de Astrofisica de Canarias, in the island of La 
Palma, as well as on data obtained with the STELLA robotic telescope 
in Tenerife, an AlP facility jointly operated by AlP and lAC. 


who detected absorption due to sodium in the atmosphere of 
ITD209458b using ITST/STIS optical observations. Today, trans¬ 
mission spectroscopy mainly comprises two observational ap¬ 
proaches; medium to high-resolution spectroscopy (5000 < R < 
100000), and low-resolution spectrophotometry {R < 200). The 
former one has a much higher sensitivity in detecting and resolv¬ 
ing very narrow spectral features, whereas the latter technique is 
able to also reveal overall ttends in the planetary spe ctrum that 
are caused b y Rayleigh scattering, for example (e.g., iPont et all 

12001120II . 

Most of the exoplanet targets for which the current genera¬ 
tion of instruments is sensitive enough to characterize the atmo¬ 
sphere are close-in gas giants. For an Earth-like extrasolar planet 
the atmospheric ring is too thin to detectably modify the host 
star’s light. The close-in gas giants, known as hot Jupiters, tend 
to have equilibrium temperature^ of 1000 K and more. Their op¬ 
tical spectra might be dominated by broad atomic alkali absorp¬ 
tion lines or, when hotter than about 1500 K, by absorption of 

* The equilibrium temperature is a theoretical temperature of the plan¬ 
etary atmosphere assuming a radiative balance between the integrated 
absorbed flux received from the star and the integrated emitted flux by 
the planet. 
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TiO according to theoretic al investigations (iHubenv et al.ll2003t 
iFortnev et al.ll200^ l2010h . Clouds of conde nsed dust particles 
are another plausible scen ario for hot Jupiters (IWoitke & Hellind 
l2Q03t iHelling et alj|2008l) in similarity to L-dwarfs whic h are 
often described by assuming a cloud layer. iFortnevI (l2005h pre¬ 
dicted that condensates even in very low densities in the probable 
atmosphere would significantly change the spectral appearance 
because of the slanted geometry under which the starlight passes 
through the planetary atmosphere. 


Optical characterizations of hot Jupiters using spectropho¬ 
tometry can already be found in the literature for multiple tar¬ 
gets. These studies show a high diversity in their optical low- 
resolution spectral appearance. 

There are robust detections of optical alkali absorp¬ 
tion features res olved in low spe ctral resoluti on R < 200 
for H D209458 b (ISing et all l2008bh . HAT-P-lb (iNikolov et al 


201 4 XO-2 dSinget al.1 120121) . and WASP-3 Ib Tsing et al 


2015 ). iMurgas et al. (2014) presented a probable detection of 
Na for WASP-43b and a tentative detection of K was given 
by INikolov et al.l (l2015h . The (low-resolution) non-detection in 
the atmospheres of HD189733b and WASP-12b can be ex¬ 
plained by scattering particles blocking the atmosphe ric heights 
that a r e possible for ato mic sodium and potassium (iPont et alJ 
l2008t ISing et aP l2013h . The non-detection of alkali absorp¬ 
tion in WASP-29b and HAT-P-32b can either be explained by 
clouds that strongly mute the spectral features or by subso lar 
abundances of so dium and potassium dGibson et al.l l20133^ . 
ISing et aP (1201 lah found significant K absorption in the atmo¬ 
sphere of XO-2b using spectrophotometry with narrow-band 
tunable filters, giving a higher spectral resolution of R ~ 800. 

A spectral slope in the optical that is indicative of 
Rayleigh scattering by a erosols has b een found for HD 189733b 


(iPont et alJl2008f 1201 3h . WASP-12b (ISing et al.ll201 3h. WA SP- 
dT l2013h . and WASP-3Tb (ISing et al.l BOll . In 


120131) .’ and WASP-3 lb ^ing et al 


6b ([Jordan et al 
contrast, the optical spectra of HD209458 b and XO-2b are 
best explained by a clear atmosphere model (ISing et al.ll2008bl 
1201 laL I 2 OI 2 I) . HD209458b also features a Rayleigh slope 
shortward of 5000 A, which is i nterpreted as scattering by 
H? mo lecules instead of aerosols dLecavelier Des Etangs et akl 

l2008bl) . 


Indications for the presence of TiO in the terminator re- 
gion was so far only found for HD29458b (iDesert et al.ll2008l) 
and WASP-12b; in the latte r case, TiH is an alte rnative solu¬ 
tion for the measurements (IStevenson et al.l |2014 . While this 
molecular absorption in HD29458b still lacks an independent 
confirmation, a follow-up HST obs ervation of WASP -12 found 
no longer any indication for TiO (ISing et alJl2013h . Nor was 
TiO was d etected in the proba bly cloud-fr ee atmospheres o f 
WASP-19b dHiiitson et al 1120131) and XO-2b dSing et al.ll2ni lah . 
Furthermore, it was r uled out for the som ewhat cooler atmo¬ 
spheres of HAT-P-lb dNikolov et al.ll2014l) . The non-detection 
of TiO in HD 189733b, HAT-P-32b, WASP-3 lb, and WASP- 


6b can be explained by 
sphere dPont et ak | 2008c [ 
INikolov rtal.ll2015h~ 


clouds or haze on top of the atmo- 
iGibson et~^l2013^ ISing et all 1201 51 : 


Here we report a search for sodium absorption in the at¬ 
mosphere of HAT-P-19b using ground-based differential spec- 
trophotometry. Th is ground-based technique was pioneered by 
iBean et ^ dioiOl) for the super-Earth GJ1214b as the search 
for chromatic transit depth variations among a set of simultane¬ 
ously observed light curves. We observed one transit event with 
the optical spectrograph OSIRIS at the Gran Telescopio de Ca- 
narias. Several studies already proved the capability of OSIRIS 


for optical transit spectrophotometry dSing et alJ l201 lai 120121 : 
IColon et aDl2012t|Murgas et al.ll2014i) . 

H AT-P-19b is a hot Ju piter with a period yery close to four 
days dHartman et al.ll20Tll hereafter HI 1). It is similar to Jupiter 
in size (1.13 Rj) and to Saturn in mass (0.29 Mj), giying it a 
yery low surface grayity of about 6 m/s^. The equilibrium tem¬ 
perature was estimated by Hll to be 1010 K. Hence, the scale 
height is ~ 600 km assuming Jupiter’s mean molecular weight 
of 2.2 g/mole. The planet transits a K-type main sequence star 
of V = 12.9 mag, M = 0.84 Mq and R = 0.82 Rq. A great 
adyantage for differential spectrophotometry is the presence of 
a near-by reference star, yery similar in brightness, just 1.5 ar- 
cminutes away and of similar spectral type. 

If present, brightness inhomogeneities on the stellar surface 
affect the deriyation of transit parameters and mimic a slope in 
the transmission spectrum. The knowledge of a long-term flux 
yariation and the stellar flux leyel at the time of the transit obser- 
yation can be used to correct for this influence. For this reason, 
we performed a two-color monitoring campaign for HAT-P-19 
using the 1.2m telescope STELLA on Tenerife oyer a time span 
of about 300 days. Similar studies can be fo und in the litera - 
ture for the hot-Jupiter host star s HD 189733 ( Pont et akl 2013). 


WASP-19 dHuitson et al l 1201 3l) . WASP- 12 dSing et al. 
HAT-P-1 dNikoloy et al.1 1201 4 WASP-6 dNikoloy et akl 
and WASP-31 dSing et alJl2015l) . 


2013), 


1201.51) . 


This paper is structured as follows: Sect. 2 describes the ob- 
seryations and the data reduction, and Sect. 3 presents the anal¬ 
ysis and results. A discussion of the results is giyen in Sect. 4, 
followed by the conclusions in Sect. 5. 


2. Observations and data reduction 

We obseryed the exoplanet host star HAT-P-19 in seryice mode 
during one transit eyent with the Gran Telescopio CANARIAS 
(GTC) located at the Obseryatorio del Roque de los Muchachos 
on the island of La Palma and operated by the Instituto de As- 
troffsica de Canarias. We monitored the star with the STELLA 
1.2m telescope of the Leibniz-Institute for Astrophysics Pots¬ 
dam (AIP) installed in the Obseryatorio del Teide on the island 
of Tenerife. 

2.1. Spectroscopic GTC OSIRIS observation 

2.1.1. Instrument setup 

During the transit obseryation we made use of GTC’s opti¬ 
cal spectrograph OSIRIS. The dispersiye element was the VPH 
R2500R, proyiding a wayelength coyerage of 5600 to 7500 A 
and a dispersion of about 1 A per binned pixel. We used the 
widest long slit ayailable at the time of the obseryation, the 5.0" 
slit, which caused the spectral resolution to be seeing limited to 
roughly R ~ 1000. We read out the two Marconi CCD detectors 
in the fastest ayailable, 500 kHz read-out full-frame mode. The 
standard 2x2 binning resulted in less than 10 seconds closed- 
shutter time between the exposures. Exposure time was 60 sec¬ 
onds. 

We chose a rotation angle of the instrument that centered (in 
dispersion direction) both the target HAT-P-19 and the reference 
star GSC 0228301197 (USNO-Bl 1246-0009792) in the slit for 
simultaneous obseryations. OSIRIS coyers its field of yiew with 
two CCD chips; we placed both objects on chip 1 to ayoid po¬ 
tential differences between chip 1 and chip 2 as error sources in 
the differential spectrophotometry. 


Article number, page 2 ofM 





























































































































M. Mallonn et al.: Transmission Spectroscopy of HAT-P-19b 


2.1.2. Observing log 

The transit observation was conducted on January 10, 2012. The 
target was observed for slightly less than five hours from 19:33 
to 00:19 UT, resulting in 238 exposures. We typically reached 
a S/N ratio of ~380 per pixel in dispersion at central wave¬ 
lengths. The transit lasted from 19:56 to 22:39 UT (first to fourth 
contact). The night was photometric with the seeing varying 
from 0.9" to 1.3". Nevertheless, the observing conditions were 
far from optimal for high-precision spectrophotometry. The ob¬ 
servation began shortly after the meridian passage of HAT-P- 
19, which occurred very close to the zenith. Unfortunately, the 
dome construction of the GTC could not be fully opened at the 
time of the observation, causing vignetting at the highest ele¬ 
vations. This vignetting was different for target and comparison 
star, clearly distracting the transit light-curve (Fig.O during pre¬ 
transit phase. Moreover, the target set quickly during the time 
series to a maximum airmass of 3.23, leading to a rapid spread 
of the point spread function with time. In Fig. [T]the evolution 
of the flux level, airmass, FWHM, and detector position relative 
to its mean is shown. The point spread function became asym¬ 
metric at about half the observing time due to a loss of focus; 
the observations were not stopped for focus adjustment to avoid 
further systematics introduced by an abrupt change of the tele¬ 
scope settings. Therefore the obtained values for the FWHM and 
pixel position need to be treated with caution because they were 
derived by fitting a Gaussian function to the spatial spectral pro¬ 
file. However, the guiding system kept the spectrum almost at 
the same position during the entire time series, it moved in spa¬ 
tial direction by less than 1 pixel. Two jumps in pixel position 
occurred, which were of subpixel value and had no effect on the 
resulting light curves. 

2.1.3. Data reduction 

The OSIRIS data were reduced by routines written in ESO- 
Midas. Standard calibration frames have been obtained at the 
same day as the science observations. Flatfield frames were ex¬ 
posed through the 5" wide long slit, the same slit as used for 
the science frames. However, calibration spectra of HgAr and 
Ne lambs were taken with a 1" wide slit to avoid unpractically 
broad emission lines. A bias subtraction of the flat, arc lamp, and 
science frames was performed using the overscan regions. The 
flatfield correction was done by a master flat created of the flux- 
weighted average of 101 single flatfield frames taken the same 
day. We noticed that the flatfield correction had almost no effect 
on the quality in terms of scatter on the final transit light-curves 
or on the absolute value of the transit parameters. For example, 
the difference in the planet-to-star radius ratio k whether a flat 
field correction was applied or not was about ten percent of the 
derived 1 cr uncertainty. 

Before the wavelength calibration we corrected for a drift of 
the spectra in dispersion direction of about 1.4 pixel over the 
length of the time series. The origin of this drift is most proba¬ 
bly a combination of telescope flexure, obvious in a shift of the 
skylines by a total of about 0.4 pixel, and a drift of the chro¬ 
matic centroids of the objects in the long slit, probably caused 
by field rotation residuals and differential atmospheric refraction 
(OSIRIS does not include an atmospheric dispersion corrector). 
For simplicity we corrected for the drift achromatically. A differ¬ 
ence in drift among the two objects has just been found for very 
high airmass, data that were excluded in the later analysis (see 
Sect. 13.11) . hence we applied the same drift correction over the 
entire spatial range. In the through-slit acquisition image, taken 


in the Sloan r band, we measured a slight mismatch of the ob¬ 
ject centroids in dispersion of ~Q.4 pixel (~0.1"), which we also 
corrected for after the wavelength calibration and spectral ex¬ 
traction. In summary, we estimate the wavelength calibration to 
be accurate to just about 0.5 A due to residual wavelength depen¬ 
dent shifts and stresses over time both in dispersion and spatial 
direction. 

In the presence of curved sky emission lines, we decided for 
a two-dimensional wavelength calibration of the whole frame 
to allow for sky estimation in spatially extended stripes distant 
to the spatial center of the object. While the choice of the sky 
stripe width minimized the light-curve scatter by up to ten per¬ 
cent depending on the wavelength region, it influenced the abso¬ 
lute transit parameters always by less than their 1 cr error values 
and did not affect the derived transmission spectrum either. The 
final width of the sky stripes spatially above and below the spec¬ 
trum was 100 pixel each. The stripes were separated from the 
object extraction stripe by a 50 pixel gap due to the very ex¬ 
tended wings of the spectral profile. The sky value at the object 
position was estimated by a linear interpolation between these 
two sky stripes, independently derived for each pixel in disper¬ 
sion direction. 

The widths of the aperture in which to extract the object flux 
were chosen to scale with the spatial FWHM of the spectral pro¬ 
file to account for the immense expansion of the point spread 
function at higher airmass (Fig.[T]). A set of scaling factors were 
tested to minimize the scatter in the white-light curve. The final 
aperture size was 7.5 times the FWHM, which ranges from 25 
to 105 pixel within the time series. Different widths of the sky 
stripes were alos tested. The scatter of the white-light curve de¬ 
creased slowly with increasing width until it reached a plateau 
at about 100 pixe ls. We determ ined whether the optimal extrac¬ 
tion technique ofE ornel (Il986l) could further minimize the dis¬ 
persion on the light curves, but as expected for very high S/N 
spectra, there was no significant difference to a simple flux sum 
in spatial direction. Instead, the simple flux sum yielded more 
robust results in the few regions with hot or dark pixels. The one¬ 
dimensional spectra were extracted in pieces of steps 10 A wide 
to account for the wavelength-dependent FWHM of the spectral 
profile. Then the pieces were stitched together in wavelength, 
forming spectra from 5617 to 7687 A. The spectra of exposure 1 
are shown in Fig.|2] 

The spectra were then divided into wavelength channels of 
a certain width (see Fig. |2] and Sect. 13.3b . For each channel the 
flux of both objects, target and comparison star, was integrated 
and used to perform differential photometry. In this way, a set of 
simultaneously observed light curves consecutive in wavelength 
was achieved. Individual photometric errors were estimated by 

( 1 ) 

with A as the area (in pixel) over which the total flux F 
(in ADU) was summed, cr the standard deviation of noise (in 
ADU) estimated in the two sky stripe areas, and g the detector 
gain of g = 1.46 e /ADU. This equation follows the photomet- 
ric error estimation of t he software package Source Extracto]0 
(iBertin & Arnoutslll996h . 


2 w w w. astromatic. net/soft ware/sex tractor 
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orbital phase 


Fig. 1. Evolution of observational parameters over the time series. From 
top to bottom: count rate of target (lower flux level) and reference star 
(higher flux level) in ADU integrated over the spectral range, airmass, 
FWHM in spatial direction at central wavelength in pixel, and spectrum 
displacement in spatial direction in pixel. 

2.2. Photometric STELLA-WIFSIP monitoring 

2.2.1. Instrument setup 

We observed HAT-P-19 with the wide-held imager WIF- 
SIP of the robotic 1.2rn twin-telescope ST ELLA on Tenerife 
jStrassmeier et al.l l2004t IWeber et al.l l2012l) . WIFSIP consists 
of a 4kx4k back-illuminated CCD with a plate scale of 0.322 
“/pixel and four read out ampliflers. It covers a held of 22 by 
22 arcminutes on the sky. The automatic scheduler of STELLA 
was set to observe HAT-P-19 on average every second night in 
blocks of three exposures in V and three exposures in I of 20 sec¬ 
onds exposure time each. We obtained 324 frames within 35 dif¬ 
ferent observing nights from December 2011 to February 2012; 
between May 2012 and October 2012, the dataset contains 1268 
frames from 82 nights. 

2.2.2. Data reduction 

The bias and flatfleld correction was made with the STELLA 
data reduction pipeline, for details see Granzer et al. (2015, in 
preparation). The following reduction steps were made with rou¬ 
tines written in ESO-MIDAS. We conducted aperture photome¬ 
try with the software package Source Extractor. Its MAG_AUTO 



wavelength [A] 

Fig. 2. Spectra of HAT-P-19 (lower flux level) and the comparison star 
(higher flux level) of the first exposure. The 50 A wavelength channels 
used for spectrophotometry are indicated by vertical dotted lines. 

option calculates an elliptic aperture individually for each image 
and object according to the second-order moments of the ob¬ 
ject’s light distribution. This method provides the flexibility to 
account for the varying observing conditions over the ten-month 
observing time. The I-band data suffered from fringing, whose 
pattern was found to depend mainly on CCD detector tempera¬ 
ture. We created master fringe maps by averaging the individual 
object-removed and smoothed science frames in groups of the 
same CCD temperature. The fringe correction was made by sub¬ 
tracting the temperature-selected master fringe map, scaled in 
amplitude to match the fringes of the individual science frames. 
The fringe residuals typically had a strength of about 20 percent 
of the original fringe pattern. 

We used the me thod of optimal w eighting of an ensemble 
of comparison stars (iBroeg et al.ll2005h to form an artificial ref¬ 
erence for differential photometry. We experimented with many 
different comparison star ensembles to verify that the observed 
photometric signal was independent of the choice of the com¬ 
parison stars. The final ensemble is formed by three stars sim¬ 
ilar in brightness and color to HAT-P-19, read out through the 
same amplifier as the target. The same ensemble was used for 
the V- and the I-band data. After averaging the three consecu¬ 
tively taken data points and subtracting a transit model using the 
transit parameter derived in Sects. ini and B.21 . the point-to-point 
scatter of the data versus a sine model (see Sect. 13.41) was about 
2.5 mmag in both filters. 


3. Analysis and results 

We fit iMandel & Agoll (l2002l) transit light-curve models to all 
transit light-curves. The main parameters of the fit are the incli¬ 
nation i, the transit midpoint Tq, the limb-darkening coefficients 
of a quadratic law ma and va, the scaled semimajor axis q/Ra 
, and the ratio of the planetary and stellar radius k - RijIRa, 
where a is the semimajor axis of the planetary orbit and R/, and 
Ra are the absolute planetary and stellar radii. The best-fitting 
parameters were determined by a least-squares Markov chain 
Monte Carlo (MCMC) approach that made use of PyAstron- 
omjQ. 


^ https://github.com/sczesla/PyAstronomy 
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3.1. Detrending and white-light curve analysis 

A white-light curve was formed by integrating the flux of the 
entire spectral range. As a very first step, we excluded the first 
20 data points that lost flux because of the obscuration of the 
primary mirror by the dome construction. Unfortunately, these 
data formed our pre-ingress measurements. The resulting light- 
curve features obviously smooth deviations from a theoretical 
symmetric transit light-curve, see Fig. [3] To model this defor¬ 
mation, we fitted a low-order polynomial to external parameters 
simultaneously to the transit model. Various polynomials to the 
external parameters time, airmass z, sky background, FWHM in 
spatial direction, pixel shift in x (dispersion direction), and pixel 
shift in y (spatial direction), and their combinations were tested. 
The pixel shift in spatial direction was measured by the fit of a 
Gaussian function to the spatial profile, the measurement of the 
pixel shift in dispersion was taken as the average of the centroid 
of six spectral features per object. 

The fit parameters of the transit model were afRA, k , i, and 
Tq with the values determined by Hll as starting values plus 
their derived errors as Gaussian priors. The eccentricity and lon¬ 
gitude of periastron were fixed to the values of HI 1. We ran ten 
chains of 1 x 10^ MCMC iterations, disgarded the initial 3x10"^ 
steps in every chain as burn-in phase, and merged all remain¬ 
ing iterations. The best-fit parameter values were derived as the 
mean of the posterior parameter distribution and the 1 cr uncer¬ 
tainties as the 68.3% of highest probability. 

To avoid overfitting of the deformations andtrends in the 
light curve, we chose the f unction that m inimized the Bayesian 
information criterion (BlC. ISchward 19781) . The BIG is similar to 
the goodness-of-fit estimation with a chi-square calculation, but 
adds a penalty term for the number of parameters in the model, 

(O — 

BIC - > -^- -H m ■ ln(n), (2) 

Z—i 0-^ 

with O in our case as the observation, M as the transit and de- 
trend model, m the number of free parameters, cr the photomet¬ 
ric error of the individual data points, and n the number of data 
points, here 175. We expect Eq. 1 to underestimate these values 
because of its ignorance of correlated noise. Therefore, we cal¬ 
culated a scaling factor to match the mean photometric error with 
the standard deviation of the photometric points. To estimate this 
factor, we only used the post-transit data (phase > 0.015) to avoid 
any bias introduced by a premature transit modeling. Then the 
BIC values for the tested low-order polynomials were calculated 
on the entire time series. Table [T] gives these values for a selec¬ 
tion of tested detrending functions together with their derived 
transit parameters, which depend significantly on the detrend¬ 
ing function choice. The final choice for the detrending function 
was a third-order polynomial with the airmass as independent 
variable because it minimized the BIC value, the scatter of the 
residuals (rms), and the amount of residual time-correlated noise 
(estimated by the y6 value, see below). 

The stellar limb-darkening was accounted for by the 
quadratic law in the transit fit. The limb-darkening coefficients 
were spherically calculated for a star of Tgjf - 5000 K and 
log g - 4.5 in high spectral resolution using PHOENIX 
stellar atmosphere models of metallicity log[M/H] = 0.0, see 
IClaret et^ (1201 3h for details. The applied values for the white- 
light curve were found by a flux-weighted average of these 
nearly monochromatic values. To analyze the white-light curve, 
we fitted for the linear limb-darkening coefficient ma and fixed 
the nonlinear term va to the theoretical value. The fitted value 
of the linear term agreed with the prediction to within about 1 



Fig. 3. Transit light-curve integrated over the entire spectral range 
(white light). The vertical dashed lines confine the data that were used 
for the analysis, see text for details. The red solid line shows the hest-fit 
model, the residuals are presented at the bottom. 


cr. This approach of fitting ma but fixing va to its theoretical 
value allows for some flexibility to account for potential sys¬ 
tematic differences between the star and its theoretical model 
spectrum used for determination of the limb-darkening parame¬ 
ters. Eurthermore, it avoids the problem of enlarged uncertainties 
when fitting for both coefficients caused by the strong correlation 
between ma and va (ISouth worth et alj|2007at ISouthworthll2008l : 

Llohnson et al.ll200^. 

We briefly comment on the systematic error introduced 
by the metallicity difference between the PHOENIX mod¬ 
els (solar metallicity) and HAT-P-19; Hll derived a value of 
log[WHL=0 .23 + 0.08. We used the limb-darkening tables of 
o are 1 (120001) and lClarel (l2004l) to estimate the systematic error 
to be on the order of 0.005 for the quadratic limb-darkening coef¬ 
ficient (the linear coefficient is a free parameter in our transit fit). 
The resulting systematic error on the derived transit parameters 
is much smaller than their uncertainties. 

The individual photometric error bars were enlarged by a 
common factor to give a reduced chi-square value of unity, 
xl - ^ -0, versus the fitted model. To account for correlated red 
noise, we derived in a second step an additional scaling factor 
as the ratio of the calculated and the theoretical standard devia- 
tion of the binned photometric residu als, often called the p factor 
(iGillon et al.l2006tlWinn et al.l2008h . We binned the residuals in 
intervals from 10 to 20 minutes (about the duration of ingress) 
in one-minute steps (original time sampling is ~1.1 minute), de¬ 
rived the factor in each case, and finally used their average as the 
P factor to enlarge the individual photometric errors. 

There are features in the noise in the post-transit light curve 
at phases later than about 0.023 visible to the eye that are be¬ 
lieved to be systematics potentially caused by differential slit 
loss. For white noise alone, the uncertainty of the radius ratio k 
should decrease approximately proportionally to the square root 
of the number of out-of-transit (OoT) data. Here it was worth try¬ 
ing whether the uncertainties actually decrease with fewer OoT 
data. And indeed, the uncertainty of the radius ratio Ak reached 
a minimum at an OoT limit of phase 0.0228. In addition, the un¬ 
certainties of the scaled semimajor axis A(fl//?A), the inclination 
A/, and the transit mid-time ATq reached their minima at about 
this number of OoT data, therefore we discarded all later phases 
from the analysis. 


Article number, page 5 oflH 





























A&A proofs: manuscript no. hatpl9_final 


Table 2. Transit fit parameter of the white-light curve. 


Parameter 

Value 

qIRa 

12.37 + 0.21 

k 

0.1390 + 0.0012 

i [deg] 

88.89 + 0.32 

Tq - 2400000 [days] 

55937.38839 + 0.00011 

Ua 

0.5736 + 0.0182 

v’A (fixed) 

0.1376 

gh [ms^2] 

5.97 + 0.61 


The results of the white-light curve analysis roughly agree 
with the work of HI 1 and are given in Tabletogether with an 
estimate of the surface gravity of the planet gh, calculated with 
Eq. 4 of I South worth et all ( l2007bl) using the value of the stellar 
velocity amplitude K from HI 1. 

3.2. Ephemeris 

We used the white-light curve of our transit measurement and 
the publicly available two complete transits from Hll, which 
were observed with the KeplerCam camera on the FLWO 1.2m 
telescope, to recalculate the ephemeris. The FLWO transits were 
reanalyzed with our algorithms to guarantee consistency in the 
analysis. A first-order polynomial over time was used for de¬ 
trending. The free transit parameter was the mid-time of the 
transit Tq. The parameters afRA, k , and i were fixed to the val- 



cycle number 


Fig. 4. Transit timing residuals versus the linear ephemeris of Eq. |3] 
The filled symbols represent the two transits of Hll and the transit of 
this work, which were used to derive the ephemeris. The open symbols 
show the amateur ETD transits. The dotted lines give the Icr uncertainty 
of the ephemeris as a function of cycle number. The dashed line shows 
the difference when the amateur data are included in the ephemeris fit, 
which yielded no significant improvement and was finally not used for 
Eq.|3] 

Table 4. Transit mid-times of HAT-P-19 and their residuals versus the 
ephemeris derived in this work. The upper part lists the values of the 
professional light curves that were used to derive the ephemeris. The 
lower part lists the values of the amateur light curves (ETD/TRESCA) 
that were tested in the analysis, but were not used in the end to derive 
the ephemeris. 


ues derived from our white-light curve, given in Table |2] The 
quadratic limb-darkening coefficients were fixed to the values 
used by Hll. The results of the mid-transit times are summa¬ 
rized in Table ID and graphically presented in Fig.|4] We experi¬ 
mented with the coefficients, using the linear coefficient as free 
fit parameter or varying both in meaningful ranges. The tests 
showed very little effect of very few seconds on Tq, which is 
well inside the errors. HI 1 used the Coordinated Universal Time 
(U TC) as time ba s e for t heir BJD calculations. As recommended 
by [Eastman et al.l (1201 Ol) . we converted all times to the Barycen- 
tric Dynamical Time (TDB) using their available online toolfl 
The newly found ephemeris is 

Tc = BJD(TDB) 2455909.326808(55) h- 4.0087844(15) A. (3) 


Mid-time Tq 

Cycle 

Residual 

Reference 

(BJD(TDB) - 24000000) 

number 

(days) 


55135.63128 + 0.00038 

-193 

-0.00011 

Hartman et al. ('2011) 

55167.70182 + 0.00040 

-185 

0.00014 

Hartman et al. (2011) 

55937.38829 + 0.00010 

7 

-0.00001 

This work 

55496.41945 + 0.00075 

-103 

-0.00255 

Naves, R 

55524.48225 + 0.00062 

-96 

-0.00124 

Muler, G. 

55528.49211 +0.00057 

-95 

-0.00017 

Ruiz, J. 

55885.27456 + 0.00050 

-6 

0.00045 

Ayiomamitis, A. 

55921.35111 +0.00075 

3 

-0.00205 

Naves, R. 

56145.84362 + 0.00051 

59 

-0.00147 

Shadic, S. 

56173.90627 + 0.00066 

66 

-0.00031 

Garlitz, J. 

56270.11642 + 0.00047 

90 

-0.00099 

Zhang, L. 


Tc is the predicted central time of a transit, N is the cycle num¬ 
ber with respect to the reference mid-time, and the numbers in 
brackets give the uncertainties of the last two digits. The refer¬ 
ence mid-time was chosen to minimize the covariance between 
reference mid-time and period. The GTC transit accordingly cor¬ 
responds to cycle 7. 

We tested whether additional transit light-curves observed by 
amateur astronomers could further improve the accuracy of the 
ephemeris. The data were made publicW available on the Exo¬ 
planet Transit Database (ETD) websit^B We selected the ETD 
transits to have uncertainties on the mid-transit time Tq of about 
one minute or less. The light curves were reanalyzed and the 
transit mid-times converted into BJD(TDB). However, we found 
that the time uncertainties of the eight ETD transits are too large 
to significantly improve the ephemeris of Eq. |3] Therefore we 
decided to use the better-documented and peer-reviewed obser¬ 
vations of Hll and our own observation. The timings of ETD 

* http://astroutils.astronomy.ohio-state.edu/time/ 

^ http://var2.astro.cz/ETD/ 


roughly agree with Eq. |2 see Fig. [H their residuals are all con¬ 
sistent with zero within less than 3 cr. 

3.3. Transit depth as a function of waveiength 

We created two different chromatic sets of light curves by bin¬ 
ning the flux of our spectra in wavelength channels of a certain 
width. At first, we chose 200 A as channel width, which yielded 
rather high S/N per photometric data point, but a poor spec¬ 
tral resolution in our resulting transmission spectrum. Then we 
binned the flux in narrower channels, which provided a higher 
spectral resolution at the cost of noisier light curves and there¬ 
fore higher uncertainties of the transit parameter. We performed 
a simple estimation of the channel width that was most sensitive 
to narrow spectral featu res with the help of a theoretical 1000 K 
transmission spectrum (iFortnev et al.l 1201^ . We compared the 
theoretically predicted transit depth in a channel centered on the 
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Na D line core with the transit depths in the adjacent channels 
for channel widths of 10 to 150 A. The function of differential 
transit depths over channel width was weighted with a simple 
square-root law to account for the S/N dependence on the width. 
For simplicity we assumed here photon noise to dominate the 
noise budget. The resulting function peaks at about 50 A. Fig¬ 
ure |5] shows the light-curve set of 200 A channel width. In the 
left column the raw light curves are shown, in the middle col¬ 
umn the detrended light curves, and in the column to the right 
the residuals versus the detrend + transit model. 

The chromatic set of light curves was analyzed simultane¬ 
ously. The fitting parameters per light curve were the four co¬ 
efficients of the airmass-dependent third-order polynomial, u\, 
and k. In addition, we fitted for a common value per set of o/Ra, 
i, and Tq with the white-light curve uncertainties of Table |2] as 
Gaussian priors. In total, the number of simultaneously fitted pa¬ 


rameters were 63 for the 10 light curves of 200 A width, and 249 
for the 41 light curves of 50 A width. The results are presented 
in Table |3] 

We verified the MCMC results of k over wavelength 
with another available transit-modeling package, JKTEBOF0 
(ISouthworth et al.ll2004l) . It uses the Levenberg-Marquardt op¬ 
timization algorithm to find the best-fitting model and in¬ 
cludes multiple options to give reliable error estimations. 
Here, we calculated the errors with a Monte Carlo simulation 
( Southworth et al.l l2005h and a residu al-permutation algorithm 
( .lenkins et al.l 2002 ; Southworthl2008h . the latter being sensitive 
to correlated noise. We adopted the higher value of both. The 
vast majority of best-fit parameters agrees within one sigma, 
and all conclusions derived in this work about the measured 
transmission spectrum of HAT-P-19b are the same for both the 
MCMC analysis and JKTEBOP analysis. On average, the uncer¬ 
tainties of JKTEBOP were larger by about 20% than the uncer¬ 
tainties of the MCMC chains. This deviation does not change 
the results of this work, and a detailed analysis of its origin is 
beyond the scope of this work. 

The transmission spectrum plotted in Pig. |6] appears very 
flat without significant outliers. We calculated the of the 
50 A data set versus a Rayleigh-scattering slope given by 
dRb/d\nA - -AH similar to the observed spectrum of 
HD 189733b (iLecavelier Des Etangs et aDl2008al) with H as the 
atmospheric scale height. We also compared the data to a ID 
solar metallicity model using a plane t-wide averaged P-T pro¬ 
file computed by iPortnev et al.l (1201 Ol) for the parameters of the 
HAT-P-19 system, see Pig. |3 Both models were binned in the 
same 50 A wavelength channels as the observed spectra. The 
cloud-free model shows a prominent sodium absorption feature 
with pressure-broadened wings, also predicted by other mod¬ 
eling wo rk for HAT-P-19b’s equi l ibrium temperature of about 
1000 K (ISeager & Sasselovll2000l: iBrownlUoOlt iHubbard et al.l 
1200 ih . A similar feature of potassium causes the upward slope 
of planetary size in the model at the red end of our wavelength 
range, the line core of potassium was not covered, unfortunately. 
The lowest of 34.7 for 40 degrees of freedom is given by a 
flat line, 36.8 was obtained for the Rayleigh-scattering spectrum 
and 48.9 was the value for the cloud-free atmosphere. 

We tried to test the influence of certain steps of our analysis 
on the obtained transmission spectrum. Por example, we shifted 
the flux channels in wavelength by fractions of their width, we 
tested different lengths of post-transit baseline and the effect of 
fixing uIRa and i to the values found by Hll. We also varied 


our treatment of the limb-darkening coefficients by fixing both 
ma and va to theoretical values, another time we fitted for both 
in the transit modeling. No test yielded a significant deviation in 
spectrum shape from the spectrum presented here in Pig.|6l 

3.4. Photometric variability of HAT-P-19 

We found a clear periodic variation of HAT-P-19 in the long¬ 
term photometry shown in PiglT] A least-squares sine fit to the 
data revealed a period of 35.5 + 2.5 days. The error value is 
drawn from the full width at half maximum of the peak in the 
periodogram in Pig. [8] The amplitude was different for the used 
filters, V and I, with 4.7+ 0.5 and 3.0+0.4 mmag provided by the 
least-squares sine fit. There are indications of varying period and 
amplitude in our data. When analyzing only the data of the ob¬ 
serving season 2012, we obtain a stronger amplitude of 6.3 + 0.6 
and 4.6 + 0.5 mmag for V and I and a slightly longer period of 
38.6 + 4.4 days. In the interpretation of the flux variation that 
is caused by starspots rotating in and out of view, these changes 
in period and amplitude are possibly related to spot evolution. 
However, the period variation is within one cr and the amplitude 
variations do not exceed 3 cr. 

The periodogram in Pig. |8]shows three additional peaks next 
to the main signal at somewhat shorter period. They loose most 
of their power or disappear entirely when subtracting the most 
significant signal (pre-whitening). Instead, in the residual data 
we find another signal above the false-alarm probability of 0.01 
at about 23 days. This period might also be related to spot evo¬ 
lution, but an in-depth period analysis is beyond scope of the 
current work. 

Hll did not mention a stellar flux variability. We reana¬ 
lyzed the HATnet photometry observed from September 2007 
to Pebruary 2008 to search for the same periodicity. The avail¬ 
able time seriefl was treated by the external parameter decor¬ 
relation technique and the trend-hltering algorithm, as is com¬ 
mon for HATnet data (see Hll and references therein). When 
we calculated a Pomb-Scargle periodogram, we found a signifi¬ 
cant period of P = 37.7+3.1 days (false-alarm probability PAP < 
0.0002), see Pig. [8] The sine fit reveals a much lower variability 
amplitude than the STEPPA data of only 1.0 + 0.3 mmag. 


3.5. Correction for starspots 


A likely reason for the periodic flux variation of HAT-P-19 are 
cool spots rotating in and out of view since the variation ampli¬ 
tude is larger in the V band than in the redder I band. Bright¬ 
ness inhom ogeneities on the stellar surface infl uence the transit 
parameters (iPont et alJl2008tlCzesla et al.ll2009ll . An occultation 
of a starspot by the transiting planet causes a bump in the light 
curve, which leads to a shallower transit and an underestimation 
of the planet size if uncorrected. Starspots unocculted from the 
planet cause the average brightness to be higher along the transit 
chord than on the rest of the stellar hemisphere, which leads to 
deeper transits and an overestimation of the planet size. Exam¬ 
ples of active host stars with frequently o bserved starspot cross¬ 
ings are HD18 9733 (see PoriLeLalJ|2()i3 a nd references therein ) 
and WASP-19 (iTregloan-Reed et al.ll2013t[Huitson et al.ll201^ . 

The white-light curve gives no evidence of a starspot cross¬ 
ing event. Possible scenarios that allow for spots along the chord 
without resolvable deformations of the transit shape such as very 
small spots or a homogeneous spot distribution along the chord 
are considered to either be negligible in their effect or unlikely. 


^ www.astro.keele.ac.uk/jkt/codes/jktebop.html 


’ http://www.hatnet.org/planets/discovery-hatlcs.html 
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Fig. 5. Chromatic set of transit light-curves of HAT-P-19b, flux integration made in 200 A wide flux channels. Left panel: the raw light curves, 
middle panel: detrended light curves, right panel: the light-curve residuals. 


We therefore assume that the transit depth is not underestimated. 
With the long-term photometry on hand, we can estimate the in¬ 
fluence of unocculted starspots on the stellar hemisphere visible 
at the moment of the transit event. This estimation depends on 
three quantities: The stellar flux level /meas at the time of the tran¬ 
sit observation, the stellar flux level without spot dimming /quiet 
, and the effective temperature of the spots Tspot. We define the 
relative variation in flux A/ by 

- /meas ~ /quiet 


with /meas as the measured stellar flux and /quiet a s the flux of 
the qu iet star without brightness inhomogeneities. iDesert et alJ 
(1201 ih described how A/ simply translates into the change in 
transit depth by 


lA -Ir 
'^meas 


'^true 


2 

true 


= aAf. 


(5) 


The proportiona lity factor a was e stimated to be -3 in the sys¬ 
tem HD 189733 (lD6sert et al.l2011|) . Other stu dies simplified Eg . 
|5]for other planetary systems by a = -1 te.g.. lBerta et al.ll201 ih . 
We cannot measure this on our own without transit measure¬ 
ments at multiple epochs, therefore we followed the simplest 
approach and set a to -1. 

A value of the flux variation A/ includes the knowledge of 
HAT-P-19’s flux level /quiet without any spots on its visible sur¬ 
face. Spots were probably also present at the time of the highest 


count rate of our light curve, therefore /quiet needed to be esti¬ 
mated. We approximated the level of a permanent flux dimming 
by the level of the va riance of the long -term light curve, an ap¬ 
proach und ertaken by Pont et al. (120131) for HD 189733 based on 
the work of lAigrain et alJ (12012 ). Thus, at the time of our tran¬ 
sit obseryation we estimate A/-0.007 in V band. Conserya- 

tiyely, we rounded this yalue to 1% flux dimming by spots. 

To estimate A/ oyer the spectral range coyered by our tran¬ 
sit obseryation, a yalue of the starspot temperature is needed. 
Howeyer, Tspot cannot easily be deriyed by photometric data be- 
cause of a degener acy between spot size and spot temperature 
(IStrassmeien 120091) . In principle, this can be broken by multi¬ 
color photometry, but our error bars on the subpercent yaria- 
tion amplitudes in V and I band cannot reasonably constrain 
the temperature contrast between spots and photosphere. The 
starspots of the exoplanet host star HD 189733, similar in stellar 
parameters to HA T-P-19, haye temperatures of between 750 K 
dSing et aTll201 Ibl) and 1000 K ( Pont et apl2008h lower than the 
spot-free photosphere. iBerdyugina l^OOSh listed seyeral early-K 
dwarfs with spots of about 1500 K temperature contrast. We cal¬ 
culated spot corrections on the transmission spectrum for three 
different temperature contrasts, 500 K, 1000 K, and 1500 K. The 
spectral energy distributions of the spots and spot-free photo¬ 
sphere were approximated by blackbody radiation. Furthermore, 
we neglected the contribution of faculae. Under these assump¬ 
tions, a flux dimming of one percent corresponds to about 2.3 
(1.4, 1.1) percent coyerage of the stellar hemisphere with spots 
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Fig. 6. Transmission spectrum of HAT-P-19b. In the upper panel the relative change in transit depth is shown for 200A wide wavelength channels, 
in the middle panel a channel width of 50 A is presented. For the lower panel two 50 A wide channels concentrated on the line center of the Na 
doublet at ~5892 A and Ha at ~6563 A. No significant variation was detected. The horizontal dotted line in the lower panel marks the average 
value, the dashed lines denote the difference corresponding to three atmospheric scale heights. 


of 500 K (1000 K, 1500 K) temperature contrast. The obtained 
correction values over wavelength range from ~1% of the 
transit depth (Ak^ ~ 0.00019) at the blue end of our spectrum 
to about 0.8% (Ak^ ~ 0.00016) at the red end. The difference in 
Ak^ according to the different spot temperatures is at the level 
of 0.1% of the transit depth at 7500 A. Our derived uncertainty 
of the transit depth k^ for the white-light curve is about 1.7%. 
Hence, the spot correction on the absolute transit depth is of mi¬ 
nor importance, and the correction on the relative changes of 
transit depth over wavelength is negligible. 


4. Discussion 

4. 1. Transmission spectrum 

We ran a simple calculation to test whether any potential vari¬ 
ations in transit depth among the narrow wavelength channels 
might be blurred by the point spread function. We assumed the 
FWHM in dispersion direction to equal the measured FWHM 
in spatial direction due to the 5” wide slit. While the FWHM 
increased to more than ten pixel toward our OoT phase limit 
of 0.0228, the average value during second and third contact of 
the transit was 4.8 pixel (1.2” on the sky). The corresponding 
light leakage weakens a potential difference in transit depth of 


the sodium 50 A channel with respect to the adjacent channels 
by ~3%. Hence, we neglect this effect here. 

Key values for an initial characterization of an exoplanet at¬ 
mosphere are the error bars of the relative transit depth measure¬ 
ments. We found the average error value to be higher than the 
scatter of the k^ values for both the 50 A and the 200 A light- 
curve sets. The reduced values of 0.9 and 0.7, respectively, 
versus a constant k^ conhrm this indicator of a slight overesti¬ 
mation of the error bars. We compared our achieved 200 A accu¬ 
racy to similar investigations present in the literature by a search 
for values of ground-based single-transit analyses using a com¬ 
parable width of the flux channels in about the same wavelength 
range. We caution that differences can occur due to different tele¬ 
scope size, target brightness, and the complexity of systematics 
in the light curves and analysis techniques. However, our average 
uncertainty of k^ for the 200 A channels of 3.6 • 10 ~"^ agrees well 
with the range of literature values 1.8 - 3.4 ■ 10 ~"^ (iGibson et al.l 
l2013albl : [Jordan et al.ll2013l;lMurgas et al.ll2014l) . One reason for 
our uncertainty to be slightly larger is certainly the lack of an 
out-of-transit baseline on one side. 

Our transmission spectrum of HAT-P-19b favors a feature¬ 
less flat spectrum, but although less likely, the Na rich, cloud-free 
atmosphere overplotted in Fig. |9] is not ruled out. If conhrmed 
by follow-up observations, a non-detection of the pressure- 
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broadened Na absorption could be explained in two ways: First, 
a cloud or haze layer does not allow probing the deeper layers 
of the atmosphere and covers the broad alkali line wings. As a 
result of the slanted viewing geometry during transmission spec¬ 
troscopy, even condens ates of rather l ow abundance can make 
the atmosphere opaque (lFortnevll2005l) . Observations of higher 
spectral resolution are needed to search for a spectroscopic fea¬ 
ture at the line core. The other explanation would be a clear at¬ 
mosphere with a depletion of atomic sodium. Sodium conden¬ 
sation into clouds of sodium sulfide, for instance, at the plan¬ 
ets cooler night side could potentially decrease the amount of 
atomi c sodium at the terminator region probed by our measure¬ 
ment (lLodderslll999t) . Another potential ca use of atomic Na de - 
pletion is ionization by stellar UV photons (iFortnev et alj|200^ . 
a scenario which is plausible because we found HAT-P-19 to be 
magnetically active. 

Among the known sample of hot-Jupiter exoplanets that have 
already been investigated by transmission spectroscopy there 
are three planets similar to HAT-P -19b in their system param- 
eters: HAT-P-lb (iBakos et alJl2007l). HAT-P-12b dHartmain et alJ 
[ 200 I . and WASP-29b ( jHellier et alj|2^ . All these planets 
have roughly the mass of Saturn, roughly the size of Jupiter, and 
an equilibrium temperature of 1000 to 1200 K. Th ey orbit solar¬ 
like ma in-sequence stars of spectral type G or K. iGibson et alJ 
(12013ah obtained an optical transmission spectrum of WASP- 
29b around the sodium line and was able t o significantly rule ou t 
a pressure-broadened absorption feature. iNikoloy et al.l (l2014l) 
obseryed HAT-P-lb in the optical regime during transit and 
detected additional absorption at the sodium core in a 30 A 
wide flux c hannel, bu t no sp ectral signature of the broadened 
line wings. iLine et aP (12013l) obseryed HAT-P-12 in the near- 
IR and found eyidence against spectral signatures of H 2 O tha 
tis theoretically expected in a clear atmosphere. All these stud¬ 
ies haye in common that they disagree with theoretical cloud- 
free atmosphere models, but could be much better explained 
by an opaque layer in the atmosphere either blocking only 
the pressure-broadened wings or the entire absorption feature 
when obseryed in low spectral resolution. Our measurements 
for HAT-P-19b agree well with these results and point toward 



54350 54400 54450 54500 

bjd +2400000 


Fig. 7. Upper panel: STELLA differential V band photometry. A simple 
sine model of the 2011 +2012 data set overplotted in green, a sine model 
of just the 2012 data in red. The vertical dotted line marks the transit 
time of our spectroscopic observation. Middle panel: The same as in 
the upper panel, here for filter 1. Lower panel: The HATnet differential 
R band photometry from 2007/2008. Black dots are the individual data 
points, red is the nightly average and green gives the best fit sine model. 


the final statement of lFortne^ (l2005h . who reported that owing 
to the slanted viewing geometry of transmission spectroscopy, 
the detection of clear atmosphere will be rare. The only ob¬ 
ject for which sodium wings, which indicate a clear atmosphere, 
have been detecte d so far is the somewhat hotter HD209458b 
(ISing et al.l2008ah . But also in the temperature regime of 2000 K 
and above, it seems to be very common to find a spectrum more 
depleted in featu res than expected by clear a tmosphere models, 
see the review bv iMadhusudhan et al.l (l2014ll . 

Future follow-up observations of HAT-P-19b with 
HST/WFC3 from 1.1 to 1.7 pm could strengthen our present 
indication for an atmosphere containing clouds or hazes by 
investigating the water absorption band at 1.4 pm. 

4.2. Long-term photometric variability 

HI 1 gathered photometry data with the HATnet telescopes from 
September 2007 to February 2008, high-resolution spectra for 
RV follow-up have been taken in August 2009 to August 2010. 
The authors did not mention any indication for stellar activity 
from photometry or spectroscopy. The RV jitter is low, 6.7 m/s, 
and might be caused by instrumental systematics. However, the 
periodic photometric variations found by us in data taken from 
December 2011 to Octobre 2012 were interpreted as starspots 
moving in and out of view because of the larger amplitude of the 
variation in B than in V. The same period of ~37 days detected in 
the HATnet data and STELLA data, consequently interpreted as 
stellar rotation period, shows an increase in amplitude from 2007 
to 2012 by about a factor of 5. Possible interpretations for this 
increase are a very homogeneous spot distribution in 2007 and a 
very inhomogeneous distribution in 2012. Another explanation 
is a change in activity level comparable with changes in the ac¬ 
tivity level during the 11-year solar cycle. Such st ellar activity 
cycles have been detected in a variety of stars, see IStrassmeiCT 
(l2009h and references therein. 

Knowing the rotation period of the star and its effective tem¬ 
perature, we have the opportunity to determine the age of HAT- 
P-19 using gyrochronology. Here, we assumed that the host star 
was not spun up by tidal interactions with its close-in planet. Ac- 



0.01 0.1 1 
frequency [cycle/day] 


Fig. 8. Lomb-Scargle periodogram for the STELLA photometry (upper 
panel) and the HATnet photometry (lower panel). Horizontal dashed 
lines denote the power value of false-alarm probability (FAP) 0.01. Ver¬ 
tical dashed lines mark the period found by the sine fit. The middle 
panel shows the STELLA periodogram after one pre-whitening step. 
Significant frequencies remain, which we attribute to the sampling pat¬ 
tern and/or spot evolution. 
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Fig. 9. Derived transm ission spectrum in c omparison to a theoretical 
spectrum calculated bv iFortnev et alJ ll2010h illustrated as a continuous 
blue line. The model values integrated over 50 A flux channel width 
are shown as blue filled circles. We also present a spectrum typical for 
Rayleigh-scattering as a red dashed line. Both models were computed 
according to the system parameters of HAT-P-19. 


cording to Eq. 15 in iHusnoo et al.l (l2012h . the timescale for syn¬ 
chronization of the rotation is above the age of the Universe and 
several orders of magnitudes higher than for systems wit h de¬ 
tected excess rotation that might be caused by tidal effec ts dPonli 
l2009t [Husnoo et al1l2012l : [Poppenhaeger &^olkll2014l) . To ap¬ 
ply the age determination from iBarnesI (1201 Oh . the effective tem¬ 
perature w as con v erted empirically to a B-V color by using the 
formula in [ReedI (Il998h . We derived B - V = 0.927 + 0.042 
from Teff = 4990 + 130 K. The errors arising from the unknown 
initial rotation period at zero age of a star between 0.12 d and 
3.4 d result in an uncertainty of 0.2 Gyr for that specific spectral 
type and rotation period. Taking all uncertainties into account 
(error in effective temperature, error in period determination, er¬ 
ror in age determination), we estimate the age of the exoplanet 
host HAT-P-19 to be 5.5^| 3 Gyr. This result refines the age esti¬ 
mation from HI 1 of 8.8 +5.2 Gyr based on stellar evolutionary 
tracks. 

Knowledge of the rotation period of the star and its spec¬ 
troscopic line broadening vsini in principle allows determin- 
ing the inclination of the star in the observer’s line of sight 
dWalkowicz & Basrill2013h . It does not provide an absolute an¬ 
gle between the orbital plane of the planet and its host star rota¬ 
tion axis, but as we know the inclination of the planetary orbital 
plane from the transit measurement to be close to 90°, any sig¬ 
nificant deviation of the stellar inclination from the same value 
indicates a misalignment of the planetary system. However, for 
very slowly rotating stars like HAT-P-19, no meaningful conclu¬ 
sions can be drawn. The vsini value of 0.7 ±0.5 km/s (taken from 
HI 1) in combination with our stellar rotation period of ~35 days 
allows for inclination angles from 10 to 90 degrees. 


5. Conclusion 

We reported GTC/OSIRIS spectroscopic observations of one 
transit of HAT-P-19b covering the wavelength range from 5600 
to 7600 A. The data allowed us to refine the transit parameters 
and the planetary orbital ephemeris. We derived an optical trans¬ 
mission spectrum of the Saturn-mass highly inflated exoplanet 
including the wavelength range of the Na D line. We extracted 
a transmission spectrum using 200 A wide wavelength channels 


(spectral resolution R ~ 30) and another using 50 A wide chan¬ 
nels {R ~ 120). The result is in both cases a flat spectrum with¬ 
out additional absorption at any specific wavelength and without 
significant trends. No sodium absorption feature of the planetary 
atmosphere was detected. 

We performed a monitoring campaign to search for a photo¬ 
metric variability of HAT-P-19 using the robotic 1.2m telescope 
STELLA in the years 2011 and 2012. We found a clear periodic 
variation with 4.7 + 0.5 mmag and 3.0 ± 0.4 mmag amplitude in 
the filters V and I and a period of ~37 days. The same periodicity 
of lower amplitude was found in the HATnet discovery photom¬ 
etry of 2007 and interpreted by us as the stellar rotation period. 
Our monitoring data covered the date of the spectroscopic tran¬ 
sit observation and could therefore be used to correct the transit 
parameters for systematics caused by spots on the visible hemi¬ 
sphere. However, the influence of the spots on the transmission 
spectrum was found to be negligible. The obtained rotation pe¬ 
riod of the host star enabled us to refine the age estimate of the 
star by gyrochronology. 

The derived transmission spectrum of HAT-P-19b favors a 
featureless gray atmosphere, but comparisons with theoretical 
models tuned to the HAT-P-19 system parameters showed that 
the achieved accuracy is not sufficient to place further constraints 
on the atmospheric composition. A model showing a Rayleigh- 
scattering slope and a solar-composition cloud-free atmosphere 
model with a pressure-broadened sodium line roughly agrees 
with the data. A transmission spectrum of higher S/N ratio can 
be obtained easily by follow-up observations of the same kind 
as presented here, since the current work is based on a single¬ 
transit observation in imperfect observing conditions. The lack 
of pre-ingress baseline made our modeling of underlying sys¬ 
tematics in the light curve (detrending) more error prone and 
hence increased our final error bars in the transmission spectrum. 
If follow-up observations confirm the flat transmission spectrum 
to higher significance, the spectrum of HAT-P-19b would share 
an important property with the spectra of other close-in gas gi¬ 
ants of_abou^JTOO_^_atmos£heric_temper^ureJnyes^ate^o 
far (iGibson et alJl2013at iNikolov et al.l 120141 : iLine et aHl2013l) : 

None of their transmission spectra can be explained by a cloud- 
or haze-free atmosphere model, they all show indications of an 
additional opacity source blocking parts of the probable atmo¬ 
sphere. 
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Table 1. Dependence of the derived transit parameters on the model function used for detrending. The notation p(f) refers to a polynomial of 
degree i of parameter], e.g., p(z^) denotes a polynomial of third degree with respect to airmass z. The model function used for the final light curve 
analysis is printed in boldface. 


model 

function 

BIC 

rms 

(mmag) 

ys 

k = Rh/RA 

qIRa 

i 

(deg) 

Piz*) 

292.1 

0.42 

1.08 

0.1406 ±0.0012 

12.32 ±0.21 

89.20 ±0.31 

Piz^) 

279.8 

0.43 

1.12 

0.1390 ± 0.0012 

12.37 ± 0.21 

88.89 ± 0.32 

Piz^) 

328.7 

0.48 

1.58 

0.1341 ±0.0020 

12.41 ± 0.43 

88.95 ± 0.61 

p{z^) + p(x') 

292.2 

0.43 

1.17 

0.1379 ±0.0012 

12.46 ± 0.25 

89.12 ±0.45 

p{z^) + p(y') 

294.4 

0.43 

1.21 

0.1382 ±0.0013 

12.23 ± 0.29 

88.95 ± 0.43 

p{z^) + p{h 

339.7 

0.48 

1.56 

0.1343 ±0.0015 

12.41 ±0.33 

89.20 ± 0.56 

Pit*) 

304.3 

0.45 

1.09 

0.1399 ±0.0014 

12.24 ± 0.27 

88.97 ± 0.46 

Pih 

319.1 

0.48 

1.36 

0.1335 ±0.0018 

12.56 ± 0.44 

89.19 ±0.64 

Pih 

316.9 

0.50 

1.51 

0.1354 ±0.0021 

12.51 ±0.36 

88.85 ± 0.59 

p(FWHM*) 

337.4 

0.47 

1.71 

0.1341 ±0.0019 

12.46 ± 0.43 

88.61 ± 0.78 

p(FWHM^) 

393.0 

0.54 

1.99 

0.1329 ±0.0027 

12.38 ±0.51 

88.95 ± 0.95 
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Table 3. Characteristics of the extracted transit light-curves. 


Wavelength 

/? 

rms 

UA 

VA 

k = VblVA 

(A) 


(mmag) 

fitted 

calculated 


5617-7616 

1.12 

0.43 

0.5736 

0.1376 

0.1390 ±0.0012 

5617-5816 

1.16 

0.73 

0.6087 

0.0988 

0.1412 ±0.0013 

5817 -6016 

1.11 

0.69 

0.6127 

0.1132 

0.1382 ±0.0015 

6017 -6216 

1.27 

0.62 

0.6233 

0.1230 

0.1390 ±0.0014 

6217 - 6416 

1.04 

0.68 

0.5964 

0.1258 

0.1392 ±0.0014 

6417 - 6616 

1.02 

0.71 

0.5663 

0.1450 

0.1392 ±0.0014 

6617 -6816 

1.15 

0.64 

0.5708 

0.1442 

0.1383 ±0.0012 

6817 -7016 

1.19 

0.68 

0.5726 

0.1495 

0.1388 ±0.0010 

7017 -7216 

1.00 

0.66 

0.5586 

0.1443 

0.1392 ±0.0012 

7217 -7416 

1.00 

0.58 

0.5155 

0.1529 

0.1407 ±0.0011 

7417-7616 

1.06 

0.59 

0.5010 

0.1570 

0.1412 ±0.0012 

5617 - 5666 

1.33 

1.21 

0.6021 

0.0989 

0.1422 ±0.0027 

5667 -5716 

1.05 

1.10 

0.5722 

0.0969 

0.1423 ±0.0023 

5717 -5766 

1.27 

1.09 

0.5748 

0.0921 

0.1408 ±0.0021 

5767 - 5816 

1.02 

0.99 

0.6060 

0.1068 

0.1380 ±0.0017 

5817 - 5866 

1.00 

0.94 

0.6293 

0.1204 

0.1358 ±0.0017 

5867 - 5916 

1.08 

1.02 

0.6034 

0.0930 

0.1394 ±0.0018 

5917 - 5966 

1.00 

1.06 

0.6470 

0.1118 

0.1359 ±0.0017 

5967 - 6016 

1.03 

1.03 

0.5886 

0.1246 

0.1379 ±0.0019 

6017 - 6066 

1.11 

1.04 

0.6163 

0.1278 

0.1377 ±0.0026 

6067 -6116 

1.22 

1.09 

0.6390 

0.1158 

0.1345 ±0.0020 

6117-6166 

1.20 

1.21 

0.6026 

0.1221 

0.1414 ±0.0022 

6167 -6216 

1.01 

1.03 

0.5902 

0.1263 

0.1389 ±0.0017 

6217 - 6266 

1.00 

0.91 

0.5626 

0.1246 

0.1397 ±0.0015 

6267 -6316 

1.18 

0.95 

0.5931 

0.1248 

0.1384 ±0.0016 

6317 -6366 

1.13 

1.05 

0.5843 

0.1307 

0.1378 ±0.0025 

6367 - 6416 

1.00 

0.99 

0.5779 

0.1232 

0.1393 ±0.0015 

6417 - 6466 

1.04 

0.90 

0.5620 

0.1225 

0.1401 ±0.0018 

6467 -6516 

1.04 

0.93 

0.5867 

0.1352 

0.1381 ±0.0019 

6517 -6566 

1.41 

0.92 

0.5425 

0.1703 

0.1402 ±0.0016 

6567 - 6616 

1.08 

1.06 

0.5467 

0.1510 

0.1379 ±0.0021 

6617 - 6666 

1.26 

1.02 

0.5639 

0.1415 

0.1391 ±0.0018 

6667 - 6716 

1.40 

0.95 

0.5864 

0.1419 

0.1364 ±0.0016 

6717 -6766 

1.31 

0.87 

0.5655 

0.1483 

0.1386 ±0.0016 

6767 - 6816 

1.06 

0.82 

0.5979 

0.1450 

0.1383 ±0.0015 

6817 - 6866 

1.00 

0.93 

0.5371 

0.1472 

0.1394 ±0.0018 

6867 - 6916 

1.43 

0.95 

0.5530 

0.1533 

0.1390 ±0.0014 

6917 - 6966 

1.00 

1.05 

0.6084 

0.1484 

0.1375 ±0.0017 

6967 - 7016 

1.15 

1.00 

0.5354 

0.1492 

0.1391 ±0.0017 

7017 - 7066 

1.05 

1.36 

0.4911 

0.1421 

0.1369 ±0.0029 

7067 -7116 

1.01 

0.95 

0.5573 

0.1498 

0.1392 ±0.0014 

7117-7166 

1.05 

1.13 

0.5190 

0.1434 

0.1393 ±0.0024 

7167 -7216 

1.02 

1.39 

0.5635 

0.1414 

0.1385 ±0.0022 

7217 -7266 

1.01 

0.97 

0.5161 

0.1494 

0.1404 ±0.0014 

7267 -7316 

1.00 

0.82 

0.4732 

0.1568 

0.1418 ±0.0016 

7317 -7366 

1.00 

0.93 

0.4958 

0.1527 

0.1403 ±0.0016 

7367 -7416 

1.04 

0.96 

0.5532 

0.1530 

0.1382 ±0.0016 

7417 - 7466 

1.00 

0.92 

0.4955 

0.1578 

0.1406 ±0.0017 

7467 -7516 

1.07 

0.81 

0.5077 

0.1572 

0.1399 ±0.0015 

7517 -7566 

1.00 

0.87 

0.4793 

0.1602 

0.1394 ±0.0015 

7567 -7616 

1.00 

0.98 

0.5008 

0.1510 

0.1432 ±0.0018 

7617 - 7666 

1.13 

1.03 

0.5131 

0.1543 

0.1394 ±0.0018 

5667 - 5866 

1.13 

0.65 

0.6002 

0.1046 

0.1395 ±0.0013 

5867 - 5916 

1.08 

1.02 

0.6034 

0.0930 

0.1394 ±0.0018 

5917 -6516 

1.25 

0.56 

0.5908 

0.1243 

0.1386 ±0.0012 

6537 - 6587 

1.31 

1.05 

0.5449 

0.1727 

0.1397 ±0.0020 

6587 - 7666 

1.14 

0.46 

0.5296 

0.1492 

0.1391 ±0.0011 




